(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 




(19) Worid Intellectual Property Organization 

International Bureau 

(43) International Publication Date (10) International Publication Number 

26 September 2002 (26.09.2002) pCT WO 02/075408 Al 



(51) InternatioDal Patent Classification^: G02B 6/34, 

H04J 14/02 

(21) International Application Number: PCT/AU02/00307 

(22) International Filing Date: 15 March 2002 (15.03.2002) 



(25) Filing Language: 

(26) Publication Language: 



English 
English 



(30) Priority Data: 
PR 3789 



16 March 2001 (16.03.2001) AU 



(71) Applicant (for all designated States except US)i RED- 
FERN OPTICAL COMPONENTS PTY LTD [AU/AU]; 
Suite 212, National Innovation Centre, Australian Technol- 
ogy Park, Eveleigh, New South Wales 1430 (AU). 

(72) Inventors; and 

(75) Inventors/Applicants ([kh- US only): STEPANOV, 



Dmitrii, Yu [RU/AU]; 7/126 Croydon Avenue, Croydon 
Park, New South Wales 2133 (AU). B.IJRYAK, Alexander 
[AU/AU] ; 53 Coronga Cresent, Killara, fJfct^XE^^Wales 
2071 (AU), 

(74) Agent: FREEHILLS CARTER SMITH BEADLE; 

Level 32, MLC Centre, 19-29 Martin Place, Sydney. NSW 
2000 (AU). 

(81) Designated States (national): AE, AG, AL, AM, AT. AU. 

AZ, BA, BB. EG, BR, BY. BZ, CA, CH. CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, Fl, GB, GD, GE, GH, 
GM, HR, HU, ID, IL, IN, IS, JP, KE. KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW. 
MX, MZ, NO. NZ. OM, PH. PL. PT, RO. RU. SD. SE, SG, 
SI, SK, SL. TJ, TM, TN, TR, TT. TZ. UA, UG, US, UZ, 
VN. YU, ZA, ZM, ZW. 

(84) Designated States (regional): ARIPO patent (GH, GM, 
KE, LS, MW, MZ. SD, SL. SZ. TZ, UG, ZM, ZW), 
Eurasian patent (AM. AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE. CH, CY, DE, DK, ES, n. FR, 

[Continued on next page] 



(54) Title: GRATING DESIGN 






< 

00 

m 
o 

(57) Abstract: A method of designing a multi-channel grating structure in a waveguide material, the method comprising the step 
of utilising a multi-channel grating design function describing an envelope of a refractive index variation defining the multi-channel 
grating structure in the waveguide material, wherein the multi -channel grating design function deviates from a periodic sampling 
^ function multiplied by a single channel grating design function. 



4.0 4.6 6.0 6.6 
I [cmj 



6.0 



wo 02/075408 Al H 



GB, GR, IE, rr, LU, MC, NL, PT, SE, TR), OAPI patent 
(BF, BJ, CF, CG, a, CM, OA, ON, GQ, GW, ML, MR. 
NE, SN. TD, TG). 

Published: 

— with international search report 



For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations ** appearing at the begin- 
ning of each regular issue of the PCT Gazette, 



wo 02/075408 





PCT/AU02/00307 



Grating Design 



Field of the invention 

The present invention relates broadly to a multi-channel grating design method and to a 
multi-channel grating structure. 

Bacicground of the invention 

Multi-channel grating stractures are typically written into photosensitive waveguides. 
The grating structure comprises a refractive index variations created in the photos^ositive 
waveguide^ which in turn determine the optical characteristics such as the reflection and 
transmission characteristics of the resulting grating structure. 

The envelope of the refractive index variation for a multi-channel grating structure is 
typically created by applying a sampling, i.e. periodic, function to a given single-channel grating 
design function. The single-channel grating profile is then typically an envelope of the refrcative 
index variations achieved by exposing the photosensitive waveguide through a suitable 
phasemask or using other interferometric techniques. 

In a vast majority of previously reported work on multi-channel gratings, a so-called 
Sinc-sampled design has been used. For the Sinc-sampling approach, an ^-channel grating 
design can be obtained by a direct in-phase summation of N identical seeding gratings [with 
k(z) — ampUtude grating ampUtude, 0 (z) - grating phase] equally spaced in the frequency space: 



Qsmc = 2m^os[(2/ -iV^-l)AAz/2] = JV^^smc[^^(Afe-2;»^)/2],sinc(Jc) s sin(x)/jc. 



and Ak is the chaimel spacing. 

This design will be referred to as "in-phase" grating design herein after by the appUcant. 

An example of tiiis design is shown in Figures 4 (c) and (d) with corresponding spectral 
characteristics as shown in Figures 4 (a) and (b). The maximum value of the refiractive index 
change required to implement this multi-channel grating design is given by a simple escpression: 



(1) 



where 



(2) 
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where An^ is the maximum refractive index change required for the single seeding grating. 
Since any photosensitive fiber used to fabricate Bragg gratings has material limits of the 
maximum achievable photoinduced refractive index change An^ this represents a linodtation on 
the maximum number of channels that can be recorded in a given fiber. Thus it is highly 
desirable to reduce a required An^, as much as possible. Also it is easy to see that in the chosen 

example (see Fig. 4) the substantial deviations from the desired (square-like in transmission; 
linear in group delay) spectral characteristics are present. It has been found that such deviations 
are always present, albeit to different degrees, if a strictly periodic sampling frmction approach 
is used« 

In another approach, one may solve an inverse scattering problem for a multi-chaimel 
grating directly (i.e. without calculating a single channel seeding profile first and then applying 
a sampling function). An example is given in Fig. 5 (a) - (d). As can be seen, the spectral 
characteristics are substanially perfect, but An^ is poorly optimised. 

At least preferred embodiments of the present invention seek to provide an altemative 
multi-channel grating design in which (i) the maximum refractive index change required as a 
fimction of &e number of chaimels is reduced when compared with the prior art grating designs 
discussed above and (ii) the resulting multi-channel gratings exhibit substantially a desired 
shape, e.g. square-like shape, in their transmission co-e£Gcient characteristics. 

Summary of the invention 

In accordance with a first aspect of the present invention there is provided a method of 
designing a multi-channel grating structure in a waveguide material, the method comprising the 
step of utilising a multi-channel grating design fimction describing an envelope of a refractive 
index variation defining the multi-channel grating structure in the waveguide material, wherein 
the multi-channel grating design fimction deviates from a periodic sampling function multiplied 
by a single channel grating design fimction. 

It has been foimd that the present invention can provide an improved grating design 
fimction when compared with prior art methods facilitating design of multi-channel gratings 
exhibiting a desired spectral characteristics (e.g. square-like shape in transmission and linear 
dq>endence in group delay). 
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In one embodiment, the method comprises calculating the grating design function, and 
the calculating comprises solving an inverse scattering problem for selected multi-channel 
spectral response characteristics, wherein each partial (single) channel response function 
describing one channel of the multi-channel spectral response characteristics includes a phase 
shift value relative to the response functions of the other chaimels. Preferably, at least one of the 
phase shift values is nonzero. In one embodiment, all of the phase shift values are nonzero. 

The method may further comprise the step of determining a set of the phase shift values 
for which an optimisation criterion is met. 

The method may comprise the step of determining a set of the phase shift values for 
which a maximum of the multi-channel grating design function amplitude is minimised. 

Altematively, the method may comprise the step of determining a set of the phase shift 
values for which a maximum diflference between a minimum and a maximum of the multi- 
channel grating design function amplitude is minimised. 

Altematively, the method may comprise the step of determining a set of the phase shift 
values for which a mean-square-deviation in the multi-chaimel grating design function 
amplitude is ntiimmised. 

The step of determining the set of phase shift values may comprise direct scanning 
through all combinations or conducting a variational analysis, or using other forms of extremum 
search numerical techniques, or a simulated annealing - Monte Carlo approach. 

In another embodiment, &e method may comprise the step of determining approximate 
values for the phase shift values. 

The det^miining of the approximate values may comprise the steps of forming a 
summation of periodic functions each describing a refiactive index variation along the 
waveguide, wherein each periodic function includes an associated phase shift value relative to 
the other periodic functions, determining a set of associated phase shift values for which an 
optimisation criteria is met, and using the set of associated phase shift values as the approximate 
values. 

Preferably, the summation of the periodic functions comprises a Fourier analysis. The 
result of the Fourier analysis may be expressed as: 
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In one embodiment, the method comprises the step of detCTnining the set of associated 
phase shift values for which a maximum of the ampUtude Q = Q{z) is minimised 

In another embodiment, tbte method comprises the step of determining the set of 
associated phase shift values for which a difference between a maximum and a minimum of the 
amplitude Q = Q{z) is minimised. 

In anottier embodiment, the method comprises the step of determioing the set of 
associated phase shift values for which a mean-square-deviation in the amplitude Q = Q{z) is 

minimised. 

The step of determining the set of associated phase shift values may comprise direct 
scanning through all combinations or conducting a variational analysis, or using other forms of 
extremum search numerical techniques, or a simulated annealing - Monte Carlo approach. 

The approximate values may be used as the phase shift values to calculate the multi- 
chaimel grating design function. 

In one embodiment, the method comprises the step of conducting a further optimisation 
process using the approximate values as seeding values for the phase shift values and using the 
results of the optimisation process for calciilating the multi-channel grating design function. 

The grating may be multi-dimensional, wherein the multi-channel grating design 
function is multi-dimensional. 

In accordance with a second aspect of the present invention, there is provided a multi- 
channel grating structure created utilising the design method defined in the first aspect of the 
present invention. 

In accordance with a third aspect of Ihe present invention there is provided a multi- 
channel grating structure, wherein an envelope of a refiractive index variation defining the multi- 
channel grating structure deviates firom a periodic sampling function multiplied by a single 
chaimel grating design function. 
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Brief description of the drawings 

Preferred forms of the present invention will now be descaibed with reference to the 
accompanying drawings. 

Figures 1 (a) -(d) show claculated spectral charcateristics and design of a multi channel 
grating structure emdbodying the present invention. 

Figures 2 (a) -(d) show claculated spectral charcateristics and design of another multi 
channel grating structure emdbodying the present invention. 

Figures 3 (a) --(d) show claculated spectral charcateristics and design of another multi 
channel grating structxire emdbodying the present invention. 

Figures 4 (a) -(d) show claculated spectral charcateristics and design of a prior art multi 
channel grating structure. 

Figures 5 (a) ^d) show claculated spectral charcateristics and design of a prior art multi 
channel grating structure. 

Figure 6 shows an experimental set up for writing a multi-channel grating structure of a 
multi-channel grating design embodying the present invention. 

Detailed description of the embodiments 

The preferred embodiment described provides a multi-channel grating design which 
exhibits desired spectral characteristics (e.g. square-like shape in transmission and linear 
dependence in group delay) and wherein the maximum refractive index change is less than 
directly proportional to the number of channels JV, thereby improving on prior art multi-channel 
grating designs. 

In the preferred embodiment, a multi-channel grating is constructed by solving the 
standard inverse scattering problem for multi-chaimel wavelength-shifted spectral 
characteristics. This has to be contrasted with prior art designs, in which a single grating design 
(seeding grating) is utilised to construct (by appling a sampling function) a multi-channel 
grating. 

Importantly, spectral response functions Hr (X) for partial single gratings are being de- 
phased with respect to* each other. In other words the inverse scattering problem may be 
presented as: 
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with nonzero relative phases 3i in the preferred embodiment. It is to be noted that 
spectral responses are being de-phased, not partial seeding gratings themselves. After solving 
the inverse scattering problem for Hr (X) the multi-channel grating design function is 
obtained and may be presented in a form: 

g(z) = KQe'^''^'^*^^ (5) 

where we explicitly retain, for illustrative purposes, the single chaimel grating design 
function The remaining factors in the expression (5) represent a "sampling" 

function, which is aperiodic. It is thus not a sampling, i.e. periodic, function anymore, but 
deviates from a sampling function. This has been found to reduce coupling of light into 
cladding modes in the resulting grating design. 

It will be appreciated by a person skilled in the art that in order to optimise the grating 
design for minimum refractive index change An^ as a function of the nxraiber of channels Ny 
suitable numerical or analytical methods can be applied. 

For relatively small number of channels N one can numerically scan through possible 
combinations of relative phases (solving an inverse scattering problem for each particular 
combination of dephasing angles i^) and selecting tiie combination which is optimal according 
to some specific selection criterion (e.g. selecting the combination which minimises maximum 
required refractive index change). 

For iV^ » 1 location of the optimal set &i is difScult Even rougih direct scanning 
through all possible sets of angles (followed by efficient numerical minimum search routines) 
quickly becomes numerically inefficient. 

To solve the optimization problem for large N we use, in an example embodiment, the 
so-called simulated annealing method - a Monte Carlo approach for minimization of multi- 
variable functions. This statistical method samples the search space in such a way that there is a 
high probability of finding an optimal or a near-optimal solution in a reasonable time. The term 
"simulated annealing" is derived from the analogy to physical process of heating and then 
slowly cooling a substance to obtain a ocystaUine structure. To start, the system state is 
initialized. A new configuration is constmcted by imposing a random displacement. If the 
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energy of the new state is lower than that of the previous one, the change is accepted and the 
system is updated. If the energy is greater, the new configuration is accepted with some 
probability. This procedure allows the system to move consistently towards lower energy states, 
yet still jump out of local n^inimal due to the probabiUstic acceptance of some upward moves. 

Another approach is to use the ^approximate equivalence between partial spectra 
dephasmg angles i9/ and partial grating relative phases (fn. Indeed, for weak gratings the JBrst 
order Bom approximation holds: 



where q{z) is a grating design function and r{fi) is a complex reflection coefficient. 
The Fourier transform (6) is a linear operation with a major property 

F(ai /^^ + fl2 f^^) = ai F /^^ + 02 F Thus, in this approximation, dephasing of partial 
gratings is equivalent to dephasing of partial spectral chaimel responses. Formally the last 
statement does not hold beyond weak grating limit. However, in practice, it is still 
approximately correct and the optimal set of angles ^ (for dephasing of partial gratings) may be 
used as a very good approximation for the optimal set of partial spectral channel angles 9i. 

Therefore, the phase shift values &i may be taken fit>m an other de-phasing grating 
design method, which is described in Patent Co-operation Treaty (PCT) patent application no. 
PCT/AU02/00160 filed on 15 February 2002, entitled 'TSlulti Channel Grating Design" assigned 
to the present applicant 

For example, in partial seeding grating dephasing described in that applicaition, a 
sampling function which periodically modulates the ampHtude of a given single-channel grating 
(seeding grating) is utiUsed, similar to prior art multi-chaimel grating designs. However, in 
addition to the periodic modulation of the amplitude of the seeding grating, different relevant 
phases (fn for each of the wavelength-seeding gratings are introduced. Accordingly, the 
resulting design function in the preferred embodiment maybe expressed as: 




(6) 



^j^'1K:oS+«+(2/-W-I)A«/2+|^^J _ 



(7) 
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where the additional phase of the grating \if-\ff{z) and the sampling amplitude Q=Q(z) 
are given by: 



g2 = 4 ^ cos(a; - )cos(w, Afe / 2 + p, ) cos(/z^ AAz / 2 + ), and 
^^^j^sina, cos(/i;Afe/2 + 



^(z) = tan 



cosa, cos(n,Afe/2 + )ff,j 



TV is even 



or 



(^-l)/2 

e^(z) = 4 ^cosa,cos(w/Afe/2 + P;) + 

/al 

4 ]^cos(a; -a^)cos(w,Afe/2 + p,)cos(«^Afe/2 + p^) + l, 



and 



\|/(z) = tan 



-1 



sm a, cos(n, AAz / 2 + P; ) 



iVis odd. 



^[^"^^'^^cosa,cos(/i,Afe/2+p, +1) 
where W; s2/-JSr-l and s2p-N-l. 

Jn the above expressions for Q(z) and t//(z) we use notations 
^ = + ^A^+i-/ ) / 2, s (^^ - ) / 2 and set (*(;y^+i)/2 = 0 for odd number of channels. Now for 
anygiven/Ztiherewillbeasetof {a,},{pj (or equivalently asetof {<|),} ) wMch miniinizes the 

ma xim um value of Q along the grating structure. By direct calculations it is straightforward to 
show that 

f;r/AJt » 
Q^dz = 2^/Ak, (8) 

for any choice of aj and P]. This expression, in turn, leads to an asymptotic formula for 
the minimum possible An^ corresponding to an "ideal" situation when Q(z)=4n and only 
phase \if{z) is nontrivially modulated: 

Anj,^yfNAn,. (9) 
We note, that, in practice, the limit Q(z) - ^fN can be reached only qjproximately. 
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Mathematically, one should solve a minimax problem and find Q^{z\a\^'^^ yfij'*''^^)fox 

which max,^^(z;a|^^'>,pj*'^'^) }=min(„^p^jmax,{Q(z;a;,p;)}. To find the optimal set for 

a relatively small number of channels one may use direct numerical scanning fiirough all 
possible combinations of the dephasing angles. 

For N»\ location of the minimizing set (a/''^'\>9/'''"^)is a nontrivial exercise. Again 
we suggest to use the simulated annealing method. 

The above described reduction strategy does not include trying to avoid 

touching the zero level at some z. However, zeros in the fibre Bragg grating (FBG) amplitude 
may lead to the increased phase errors (appearance of phase jumps) and shoxild be avoided. 
Thus, arguably, a better minimization strategy is minimizing its maximum deviations of Q{z) 

along z firom the theoretical limit level of 4n . Mathematically this may be formulated as 

This approach may be implemented by using the same simulated annealing algorithm 
described above. 

Another approach will now be described, in which optimization by the fimctional 
minitnisation (variational approach) is utilised. The key property of this embodiment is that it 
relies on estimate of some integral functional rather than time-consuniing numerical scanning in 
z. Quantitatively, proximity of Q{z) to the theoretical Umit 4n can be characterised by mean- 
square-deviation, 

^Q = ^{{Qiz)-Qy)> (10) 

where (/(z)) s / = f(z)dz Ideal optimisation of Q(z) corresponds to the 

achievement of the average and the zero mean-square-deviation fix)m this average 

value. Using expression (4) and assuming Q « 4n , one can show that 

AQ^^i2-E)EN, (11) 
where 
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cos(a„ "•a„^,,.^,)cos(a; -a^ -p^ +p„ -p,^,_^|) . 



For finding noLmima of t^Q the most efficient strategy is again the use of the simulated 
annealing method. 

The advantage of the optimisation based on the functional noinimisation compared with 
direct scanning is the speed: integration over z is carried out analytically which saves lots of 
computer time. For odd number of channels variational optimisation leads to a sampling 
function without zeros in amplitude (similar to maximum deviations minimisation approach). 
For even number of channels variational optimisation leads to a sampling function with zeros in 
amplitude (similar to maximum minimisation approach). 

The calculated angles (a/''^^\/?/'^^) maybe used directly for partial spectra dephasing of 
multi-chaimel grating designs to derive a multi-chaimel grating design function embodying the 
present invention, which deviates fix)m a periodic sampling function multipUed by a single 
channel grating design function. In an embodiment of the present invention, phase shift values 
(see equation (7)) corresponding to maximum minimisation ^proach for a four-channel 
grating design were determined utilising the above optinoisation calciilations, in which ay = 
0.575971:, - 0, and fii- fii-O was chosen. 

The phase shift values thus determined were lhaa utilised as phase shift values &i for 
the inverse scattering problem as part of this embodiment of the present invention (see equation 
(4)). The envelope of the refiractive index variation and transmission characteristics are shown 
in Figures 1(b) and 1(a) respectively. It can be seen that the grating design of the preferred 
embodiment retains all desirable spectral characteristics with an almost absolute accuracy. In 
Figure 1(d) the envelope of the refiractive mdex variation of Figure 1(b) has been normalised to 
the single channel grating design function used, over a limited range, to illustrate the deviation 
firom a periodic sampling function used in prior art designs. The calculated time-delay 
characteristics of this embodiment is shown in Figure 1(c). 

Another embodiment of the present invention using partial spectra dephasing of multi- 
channel grating design (based on difference minimisation approach) is shown in Figures 2 (a) - 
(d). The envelope of the refiractive index variation and transmission characteristics are shown in 
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Figures 2(b) and 2(a) respectively. It can be seen that the grating design of this embodiment also 
retains all desirable spectral characteristics with an ahnost absolute accuracy. In Figure 2(d) the 
envelope of the refractive index variation of Figure 2 (b) has been normalised to the single 
channel grating design function used, over a limited range, to again illustrate the deviation from 
a periodic sampling frmction used in prior art designs. The calculated time-delay characteristics 
of this embodiment is shown in Figure 2(c). All said about the grating design shown in Figure 1 
stays valid, but in addition zeros ui grating ^odisation profile are avoided. 

Both of the designs embodying the present invention provide about 40% reduction in 
maximum An^ in comparison with Sinc-sampling or in-phase inverse scattering-based prior art 
designs (See Figs. 4 or 5) and at the same time have substantially ideal characteristics, similar to 
in-phase inverse scattering-based prior art designs (see Fig. 5). 

It will be appreciated by a person skilled in art that the phase shift values determined 
utilising the "partial seeding grating dephasing" in the embodiments described above may, in 
other embodiments, be utilised only as an initial set of phase shift values for solving the inverse 
scattering problem for the multi-chaimel wavelength shifted spectral responses. The method of 
such embodiments then fiirther comprises the step of finetuning of dephasing angles usrug of 
the initial set of phase shift values as seeding values. Ih that way, a grating design is fiirther 
optimised by fiirther reducing the maximum refiractive index of a given grating design (with 
seeding dephasing values). The standard nimierical minimum search techniques (in combination 
with efficient inverse scattering numerical methods) are utilised for this purpose. 

An example of such a finetuned optimisation embodiment is shown in Figure 3. The 
envelope of the refractive index variation and transmission characteristics are shown in Figures 
3(b) and 3(a) respectively. It can be seen that the grating design of this embodiment again 
retaixis all desirable spectral characteristics with an almost absolute accuracy. In Figure 3(d) the 
envelope of the refractive index variation has been normalised to the single channel grating 
design fimction used, over a limited range, to illustrate the deviation from a periodic sampling 
fimction used in prior art designs. The calculated time-delay characteristics of this embodiment 
is shown in Figure 3(c). One can see that the maximum refiractive index change of this design is 
sligjitly (5%) smaller, than the maximum refractive index change of the design shown in Figure 
1. Spectral characteristics still remain as good as for the design of Fig. 1. 
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As can be seen from Figures 1 to 3, the implementation of the mnlti-channel grating 
design of the preferred embodiment in a grating structure requires grating writing apparatus 
with high spatial resolution to be utilised. Therefore, in a grating writing apparatus relying on 
photo induced refractive index changes, the apparatus preferably comprises a beam focusing 
means to reduce the size of the beam in the core of the photosensitive waveguide. 

Figure 6 shows an example experimental set up 50 for writing a multi-channel grating 
52 into an optical fibre 54. The experimental set up 50 comprises an interferometer 56 which 
includes a first acousto-optic modulation 58 being operated under an acousto-optic wave of a 
first frequmcy Qi, as indicated by arrow 14. An incoming light beam 60 is incident on the first 
acousto-optic modulator 58 under a first order Bragg angle. The operating conditions of the 
acousto-optic modulator 58 are chosen such that the modulator 58 is under driven, whereby 
approximately 50% of the incoming beam 60 is diffracted into a first order beam 62, and 50% 
passing through the acousto-optic modulator 58 as un-diffracted beam 64. The un-diflBracted 
beam 64 is incident on a second acousto-optic modulator 66 of the interferometer 56 under a 
first order Bragg angle, whereas the beam 62 is not. Accordingly, the beam 62 passing through 
the second acousto-optic modulator 66 without any significant loss. 

The second acousto-optic modulator 66 is operated under an acoiisto-optic wave of a 
frequency Q2» which propagates in a direction opposed the direction of the acousto-optic wave 
in the first modulator 58 as indicated by arrow 68. After the second acousto-optic modulator 66 
the first order diffracted beam 70 and the beam 62 are frequency shifted in the same direction 
(e,g. higjier frequency), but by different amounts i.e. Qi v Qa- 

The beams 62, 70 are then brou^t to interference utilising an optical lens 72, and the 
resulting interference pattern (at numeral 74) induces refiractive index changes in the 
photosensitive optical fibre 54, whereby a refractive index profile, i.e. grating structure 52, is 
created in the optical fibre 54. 

In Figure 6, the optical fibre 54 is translated along the interferometer at a speed v, as 
indicated by arrow 74. 

It will be appreciated by a person skilled in the art that the experimental set iip 50 shown 
in Figure 6 can be utilised to write a multi-channel grating structure of a multi-channel grating 
design embodying the present invention through suitable control of the first and second acousto- 
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optic modulators 58, 66, in conjunction with a smtable control of the speed v at which the 
optical fibre 54 is translated along the interferometer 56 at any particular time. The high spatial 
resolution required to implement the multi-channel design of the preferred embodiment is 
achieved in the set \xp shown in Figure 6 by utilising optical lens 72, with the practical limit of 
the beam size in the focal plane preferably being of the order of the waveguide core size. 

It will be appreciated by a person skilled in the art that numerous variations and/or 
modifications may be made to the present invention as shown in the specific embodiments 
without departing fi?om the spirit of scope of the invention as broadly described. The present 
embodiments are, therefore, to be considered in all respects to be illustrative and not restrictive. 

For example, multi-channel gratings can be created on the basis of the multi-channel 
grating design of the preset invention using various known grating creation techniques, 
including one or more of the groiq) of photo-induced refractive index variation in photo 
sensitive waveguide materials, etching techniques including etching techniques utilising a 
phasemask, and epitaxial techniques. 

In the claims that follow and in the summary of the invention, except where the context 
requires otherwise due to express language or necessary implication the word "comprising" is 
used in the sense of "including", i.e. the features specified may be associated with fiirther 
features in various embodiments of the invention. 
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Claims 

1 . A mefhod of designing a multi-channel grating structure in a waveguide material, 
the method comprising the step of: 

- utilising a multi-channel grating design fimction describing an envelope of a refractive 
index variation defining the multi-channel grating structure in the waveguide material, 

wherein the multi-channel grating design function deviates from a periodic sampling 
function multiplied by a single channel grating design function. 

2. A method as claimed in claim 1, wherein the method comprises calculating the 
grating design function, and the calculating comprises: 

- solving an inv^e scattering problem for selected multi-channel spectral response 
characteristics, 

wherein each partial (single) channel response function describing one channel of the 
multi-channel spectral response characteristics includes a phase shift value relative to the 
response functions of the other channels. 

3. A method as claimed in claim 2, wherein at least one of the phase shift values is 
nonzero. 

4. A method as claimed in claim 3, wherein all of the phase shift values are 
nonzero. 

5. A method as claimed in any one of claim 2 to 4, wherein the method further 
comprises the step of determining a set of the phase shift values for which an optimisation 
criterion is met. 

6. A method as claimed in claim 5, wherein the method comprises the step of 
determining a set of the phase shift values for which a maximum of the multi-channel grating 
design function ampUtude is minimised. 

7. A method as claimed in claim 5, wherein the method comprises the step of 
determining a set of the phase shift values for which a maximum difference between a 
minimum and a maximum of the multi-channel grating design function amplitude is nunimised. 
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8. A method as claimed in claim 5, wherein the method comprises the st&p of 
deteraiining a set of the phase shift values for which a mean-square-deviation in the multi- 
channel grating design function ampUtude is minimised. 

9* A method as claimed in any one of claims 5 to 8, wh^eia the step of determining 
the set of phase shift values comprises direct scaiming through all combinations or conducting a 
variational analysis, or using other forms of extremimi search numerical techniques, or a 
simulated annealing - Monte Carlo approach. 

10. A method as claimed in any one of claims 1 to 4, wherein the method comprises 
tiie step of determining ^proximate values for the phase shift values. 

11. A method as claimed in claim 10, wherein the determining of the approximate 
values comprises the steps of: 

- forming a summation of periodic functions each describing a refractive index variation 
along the waveguide, wherein each periodic ftmction includes an associated phase shift value 
relative to the other periodic ftmctions, 

- detemiining a set of associated phase shift values for which an optimisation criteria is 
met, and 

- using the set of associated phase shift values as the approximate values. 

12. A method as claimed in claim 11, wherein the summation of the periodic 
ftmctions comprises a Fourier analysis. ^ 

13. A method as claimed in claim 12, wherein the result of the Fourier analysis is 



14. A method as claimed in claim 13, wherein the method comprises the step of 
determining the set of associated phase shift values for which a maximum of the amplitude 
Q = Q(z) is minimised. 




IS. A method as claimed in claim 13, wherein the method comprises the step of 
detemiining the set of associated phase shift values for which a difference between a maximum 
and a minimum of the amplitude Q = Q(z) is miiiimised. 
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16. A method as claimed in claim 13, wherein the method comprises the step of 
detemiining the set of associated phase shift values for which a mean-square-deviation in the 
amplitude Q = Q(z) is minimised. 

17. A method as claimed in any one of claims 13 to 16, wherein the step of 
detemiining the set of associated phase shift values comprises direct scanning through all 
combinations or conducting a variational analysis, or using other forms of extremum search 
numerical techniques, or a simiilated annealing - Monte Carlo approach. 

18. A method as claimed in any one of claims 12 to 17, wherein the approximate 
values are used as the phase shift values to calculate the multi-channel grating design fimction. 

19. A method as claimed in any one of claims 12 to 17, wherein the method 
comprises the step of conducting a ftirttier optimisation process using the approximate values as 
seeding values for the phase shift values and using tihe results of the optimisation process for 
calculating the multi-channel grating design ftmction. 

20. A method as claimed in any one of the preceding claims, wherein the grating is 
multi-dimensional, and wherein the multi-channel grating design ftinction is multi-dimensional. 

21. A multi-channel grating stmcture created utilising the design method as claimed 
in any one of the preceding claims. 

22. A multi-channel grating stmcture, wherein an envelope of a refractive index 
variation defining the multi-channel grating stmcture deviates bom a periodic sampling 
ftmction multiplied by a single channel grating design ftmction. 
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